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1. Introduction 

With the restart of the LHC at an unprecedented 13 TeV center-of-mass energy, the search 
for beyond-standard-model particles continues. For a successful achievement of the goals of this 
second run, precision is vital from both theory calculations and experimental measurements. Due 
to the inevitable hadronic environment involved at the LHC, where for example most Higgs studies 
are overwhelmed by QCD background, analytic estimates of QCD observable cross-sections will 
continue to play a central role leading either to the systematic improvement/tuning of Monte Carlo 
event generators, or to the development of better methods of background elimination. 

The resummation of large logs, resulting from the reaVvirtual mis-cancellation of soft and/or 
collinear singularities in the matrix element, is perhaps the most challenging QCD perturbative 
aspect when one attempts to make an estimate of the cross-section of a given observable V. For 
observables that are sensitive to emissions in the entire angular phase space, and which are termed 
“global” observables, the resummed distribution maybe cast into the general form: 

a(F) °c exp (Lgi(a,L) +g 2 {a,L) -f a,g 3 {a,L) + afg^ia^L) + ■■■), (1.1) 

where L is the large log of the observable V. The functions g\, g 2 , ■■■, respectively resum leading 
logs (LL), next-to-leading logs (NLL),_ 

While the development of the resummation programme for global observables has seen sub¬ 
stantial progress in recent years, achieving up to NNNLL accuracy (i.e. up to gn in eq. (0), 
e.g. for the C-parameter distribution ||T|], and even semi-automatic resummation to NLL |Q] and 
NNLL [Q], progress in the resummation of “non-global” observables [^, ^ has been very limited. 
Non-global observables are those which are sensitive to emissions in restricted regions of the an¬ 
gular phase space, and as a result their distributions suffer from non-global logs (NGLs). There 
are several important observables that are non-global and that are widely used in studies relevant 
to new physics searches, such as jet mass and single hemisphere observables. 

There are two main reasons that have long jeopardised progress in the resummation of NGLs. 
Firstly, the treatment of cascade gluon branching is extremely cumbersome within perturbation 
theory. Secondly, the phase-space integrations that one has to perform are prohibitive at higher 
orders due to non-iterative geometry involved in the calculation. A practical solution to the above 
hindrances is resorting to the large-Nc approximation, with Nc being the number of quark colours, 
which amounts to discarding non-planar Feynman diagrams in the calculation of amplitudes 
of soft gluon emissions, leading to considerable simplifications. Specifically fhis approximafion is 
equivalenf fo fhe leading-order expansion of fhe colour faclor Cp = Nc/2 — l/2Nc ~ Nc/2. Fur- 
fhermore, a convenienf approach fo deal wifh fhe mulfi-dimensional infegrafions is restoring to 
numerical Monfe Carlo mefhods. The numerical resummafion of NGLs af large Nc was firsl per¬ 
formed in fhe pioneering work of Dasgupfa and Salam 0. 

There has recenfly been an increasing inferesf in liferafure in fhe calculation of NGLs bofh 
af fixed order and fo all orders. In ref. [^, Rubin numerically evaluated NGLs af large Nc for 
bofh filtered Higgs jef mass and inferjef energy flow up fo sixfh order in a^. Haifa and Ueda [^] 
performed a numerical resummafion of NGLs based on fhe Weigerf equalion 0, accounting for 
NGLs af finite Nc fo all orders. A similar equalion lhal was developed by Banfi, Marchesini and 
Smye, fhe BMS equalion [|T^], and whose solution accounfs for fhe all-orders resummafion of 
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NGLs at large Nc, was the subject of study by Schwartz and Zhu in ref. O- where the analytic 
solution up to fifth order was achieved. 

The aim of this work is to address the question of how the accuracy of resummation of NGLs 
is affected by the large-Nc approximation. For this, and other reasons, we perform the analytic 
calculation of NGLs at finite Nc up to fifth order. In the next section we define fhe observable fhaf 
we use fo illusfrale fhe calculafion of NGLs, namely fhe single-hemisphere mass disfribufion in 
—>• qq. In secfion 3 we show fhe resulfs for fhe NGLs up fo fiffh order, which are fhen used fo 
make an ansfaz for fhe all-orders (partial) resummation of NGLs into an exponential form. We also 
compare, in the same section, our findings fo fhose reporfed af large Nc by Schwarfz and Zhu [jn]]. 
In secfion 4 we perform a comparison wifh fhe numerical resulfs obfained by Dasgupfa and Salam 
0] and discuss fhe implicafions of our resulfs. Finally we summarise and give fulure directions of 
Ibis work in secfion 5. 

2. Observable and kinematics 

We are inferesfed in fhe calculafion of NGLs af finife Nc up fo fiffh order. For illusfrafive 
purposes we choose fo sfudy fhe simple process qq accompanied by fhe emission of sofl 

energy-ordered gluons ki, as depicfed in figure |i[ The quark and anfi-quark direcfions defermine 



Figure 1: Diagrams for gluon emission from an outgoing qq dipole, relevant for NGLs calculation. 


two back-to-back hemispheres and We consider for measurement the right hemisphere 
Jf^R and calculate its normalised invariant mass p defined by: 

P = (p<i+Y,k\ ( 2 . 1 ) 

\ / i^-^R 


where Q is fhe cenfer-of-mass energy and Pq and ki are fhe four-momenfa of fhe quark and gluons, 
respectively. 

The infegrafed hemisphere-mass disfribufion, normalised fo fhe Born cross-section rto, is: 

r ^^^P'= '^ + <^iiP) + <^ 2 ip) + ---, (2.2a) 

Jo rtodp' 

(ym=Y, [ > (2-2b) 

X 2 


where = W^{ki,k 2 ,...,km) represenfs fhe eikonal squared amplifude, normalised fo fhe 

Bom squared amplifude, for fhe emission of m energy-ordered gluons and d(p,n is fhe corresponding 
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phase space. The sum over X accounts for all possible real/virtual configurations of the radiated 
gluons. The measurement operator plays the role of an event selector, i.e. it forbids real 
emissions into M’r which contribute more than p to the hemisphere mass. 


In order to be able to compute the cross-section ( |2.2[ ), there are several issues that need to be 
addressed. First, one must evaluate the gluon-emission squared amplitudes including all 

the possible real-virtual configurations X at each order in the perturbation series (our aim is up to 
fifth order). The calculation of such amplitudes is non-trivial due to the complexity of the colour 
algebra involved as well as the factorially growing number of Feynman diagrams that one has to 
account for at each order. Further details about the computation of these eikonal amplitudes, which 
involves using the Mathematica package ColorMath [|I^, are to be found in our work in refs. 

00 - 

The second task is to apply the measurement operator according to the various real-virtual 
gluon configurations in order to extract the appropriate phase-space region of integration for each 
gluon. Doing so the final task is to perform the relevant multi-dimensional integrations in eq. 

At fourth order, for instance, the integral is seven-dimensional and was performed semi- 
analytically. 


3. Non-global logs up to fifth order at finite Nc 

We can write the integrated hemisphere mass distribution as: 

a(p) = a^(p)xa^°(p), 

a^(p) =exp(-CFaiT2), 


(3.1a) 

(3.1b) 


where L = ln(l/p) and as = a^/TT. We have factorised the distribution into the product of a 
Sudakov form factor a^, that resums double logs originating from soft-collinear primary emissions, 
times a non-global factor that resums single logs originating from soft wide-angle secondary 
correlated emissions. 

We express the non-global factor as a series in the coupling starting from second order, where 
NGLs first appear, up to fifth order as follows: 


=1 - —CfCaC2 + —CfC^Cs - ^ f yCFCAC4 - yCpC^G 


i2/^2 


^c2ciC2C3 + ^CFCiC5 


O / Cf 1 


+ 


'i J ) 


(3.2) 


with L = cCgL, Ca = Nc, and a and p are constants that are yet to be determined. 

To compare our result (^) with that obtained in ref. [ [n[ ] at large Nc, by means of analytic 
solution to the BMS equation, we simply make the substitution Cp —> Nc/2, leading to: 


^''°(p)=1-^(NcI)' + ^(NcL)3 + 


34560 


.4 _ / -2 

(NcL)V - 




288 

(3.3) 


which exactly agrees with the result arrived at in ref. | pl] ] up to fourth order. Furthermore, at fifth 
order we can extract the value of the undetermined constant a by comparison with the result in ref. 
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[ pl| ] and we obtain a = 17/2 + (^ 2 ^ 3 /Cs- We can also further make an anstaz for the constant p 
based on the pattern of zeta functions observed at previous orders: p =21^2^31^5- 

In order to check the impact of finite-Nc corrections (relative to large-Nc result) on the distri¬ 
bution, we compare the result at large Nc (eq. to that at finite Nc (eq. (^)). We find that at 

fourth order the size of finite-Nc result constitutes merely an ^(1.5%) correction to the large-Nc re¬ 
sult. This observation is in accordance with that made in ref. ||^ by means of numerical evaluation 
of NGLs at finite Nc through a solution to the Weigert equation. 

We note that the series of NGLs in eq. adheres a pattern of an expansion of an exponential 
function: 


a^^(p) =exp ( -^CfCa ?2 + ^CpCiCa - 


29 /Cf 1 
^"2 


+ 


+ fcFClC5 




= exp 


r 2 r3 tA 

-^CfCaC2 + ^CfCaCs - 


-h 


5! 


17 


CFClC5+2C2ciC2?3 




25 


+ 1 


a" 


CFCiC4 + ^c2ciCl 


-h 


a" 


(3.4) 


where we substituted the values of a and p based on the observations discussed above. The first 
form of the above exponential explicitly displays the finite-Nc corrections, whilst the second form 
focuses on disclosing the pattern of the Zeta functions. 

To determine the convergence of the series in the exponent of eq. (Q, we plot in figure | 
fhe ratio a'^^/exp(a|^^), wifh = —L^/2! x CFCAC 2 ^ being fhe second-order NGLs function 



Figure 2: Plot of the ratio (7^*^/exp((J2*‘^). 


in the exponent. In figure ^ we show fruncafions of fhe series af fhird, fourth and fiffh orders. If is 
clear from fhe curves fhaf fhe fhird, fourfh and fiffh-order ferms in fhe exponenf form a significanf 
1^(30%) confribufion (parficularly af larger values of L), meaning fhaf fhe series converges very 
slowly. If also means fhaf more terms of fhe series in fhe exponenf are needed for a phenomenolog¬ 
ically reliable esfimafe of fhe all-orders behaviour of fhe disfribufion. 
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4. Comparison to numerical results at large Nc 


In this section we compare the results we obtained for the resummed NGLs, the exponential 
form with the numerical all-orders resummed result obtained by Dasgupta and Salam at large 
Nc via a Monte Carlo approach 0]. Their numerical result is parameterised as follows 0: 


=exp 


-CfCa 


3 


l + (0.85CAt)2 
1 + (0.86CaO''^^^ 


(4.1) 


with the evolution parameter t given by: 


-In 


dTTjSo 1 - 2j8oa,L ’ 


(4.2) 


and j8o is the leading-order coefficient of the QCD jS function. We show in figure |3| (on the left) 
a plot of the functions (eq. ([4.1[)) and (eq. (3.4)) with various truncations of the NGLs 
series in the exponent of eq. 0. We also show in the same figure (on fhe righf) a plof of fhe ratio 
of fhe fwo funcfions. 




Figure 3: Left: Plot of the functions ffp® and with various truncations of the NGLs series. Right: Plot 
of the ratio 


It is clear from the plots, particularly the right-hand-side one, that as one adds more terms 
in the exponent of eq. (3.4), one obtains larger intervals (starting from t = 0 and spanning over 
large values of t) over which there is an agreement between the analytical form (^^) and the 
parameterised form of the Monte Carlo output (|4. l| ). This observation means that it may suffice fo 
compufe jusf a few more higher-order ferms in order fo obfain agreemenf befween fhe fwo funcfions 
for a phenomenologically sufficienf interval of t. If is worfh noting fhaf fhe second-order resulf has 
fhe peculiar fealure fhaf if represenfs fhe besf fif fo fhe all-orders resulf over fhe full range of t 
considered. 


5. Summary and outlook 

In fhis work we addressed fhe calculation of NGLs af finife Nc af single log accuracy up fo 
fiffh order for hemisphere mass disfribufion in 2 jefs. This was achieved fhrough a brute- 

force mefhod in which we infegrafed eikonal squared amplifudes over fhe appropriafe phase space. 
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We observed that the obtained series of NGLs suggest a possible resummation into an exponential 
form. When expanded to leading order in colour, our results exactly reproduce those obtained at 
large Nc by means of the analytical solution to the BMS equation. The results we obtained agree 
with the statement made in ref. [||] that finite-Nc contribution forms a small correction to the large- 
Nc result, meaning that the large-Nc approximation is a good one, at least in the context of 
collisions. 

Our next task in this work is to go beyond fifth order in the calculation of the series of NGLs 
in an attempt to confirm fhe sfrucfure of fhe resummation into an exponential form. In addition we 
plan to extend this work by analytically investigating the effect of jet clustering on the hemisphere 
mass as well as other jet-shape distributions. 
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